arXiv:1501.07897v2 [astro-ph.IM] 19 Feb 2015 


The Whole is Greater than the Sum of the Parts: 
Optimizing the Joint Science Return from LSST, Euclid and 

WFIRST 


August 8, 2016 


B. Jain/ D. Spergel/ R. Bean, A. Connolly, I. Dell’antonio, J. Frieman, E. Gawiser, N. Gehrels, 
L. Gladney, K. Heitmann, G. Helou, C. Hirata, S. Ho, Z. Ivezic, M. Jarvis, S. Kahn, J. Kalirai, 
A. Kim, R. Lupton, R. Mandelbaum, R Marshall, J. A. Newman, S. Perlmutter, M. Postman, 
J. Rhodes, M. A. Strauss, J. A. Tyson, L. Walkowicz, W. M. Wood-Vasey 


Contents 

1 Where will we be in 2024? |2] 

2 Background on missions 

2.1 WFIRST. 

2.2 LSST. 

2.3 Euclid. 

3 Enhanced science 

3.1 Weaklensing. 

3.2 Large-scale structure. 

3.3 Galaxy clusters. 

3.4 Supernovae. 

3.5 Strong tensing. 

3.6 Photometric redshifts. 

3.6.1 The impact of Euclid and WFIRST near-infrared data on LSST photometric redshifts 

3.6.2 Mitigating systematics with WFIRST and Euclid spectroscopy. 

3.7 Stellar science. 

3.8 Galaxy and quasar science. 

3.9 The variable universe. 

4 Coordinated observations 

4.1 Measuring dark energy with strong lens time delays. 

4.2 Probing dark matter with grism spectroscopy of lensed AGN. 

4.3 Deep drilling fields. 

4.4 Coordinated supemovae observations . 

’ bj ain @ physics.upenn.edu 
^dns@ astro.princeton.edu 


1 






















4.5 Variable universe: follow-up on novel targets . 

5 Data management and analysis challenges |24] 

5.1 Multi-resolution photometry. [24] 

5.2 Joint simulations . [25] 

5.2.1 Cosmological simulations. [^ 

5.2.2 Instrument and pipeline simulations. [2^ 

6 Conclusion [27] 


1 Where will we be in 2024? 

Astronomy in 2024 should be very exciting! LSST and Euclid, which should each be in the midst of their 
deep surveys of the sky, will be joined by WFIRST. With higher resolution and sensitivities than previous 
astronomical survey instruments, they will reveal new insights into areas ranging from exoplanets to the 
nature of dark energy. At the same time, JWST will be staring deeper into the early universe than ever before. 
Advanced LIGO should be detecting frequent collisions between neutron stars. ALMA will be operating at 
all of its planned frequencies, and the new generation of very large optical ground based telescopes should be 
revolutionizing ground-based optical astronomy. In parallel, advances in computational capabilities should 
enable observers to better exploit these complex data sets and theorists to make detailed time-dependent 
three-dimensional models that can capture much of the physics needed to explain the new observations. 

The focus of this report is an exploration of some of the opportunities enabled by the combination of 
LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade’s astronomy. 
The sum of these surveys has the potential to be significantly greater than the contributions of the individual 
parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high- 
resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar 
and galactic data have the potential of yielding new insights into topics ranging from the formation history of 
the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this 
multi-instrument data set is a challenging technical task: for much of the science, we will need to combine 
the photometry across multiple wavelengths with varying spectral and spatial resolution. Coordination will 
be needed between the LSST, Euclid, and WFIRST projects in order to understand the trades between 
overlapping areal coverage, filter design, depth and cadence of the observations, and performance of the 
image analysis algorithms. We will need to provide these data to the community in a highly usable format. 
If we do not prepare the missions for this task in advance, we will limit their scientific refum and increase 
fhe cosf of fhe evenfual efforf of fully exploiting fhese dafa sefs. The goal of fhis reporf is fo idenlify some 
of fhe science enabled by fhe combined surveys and fhe key technical challenges in achieving fhe synergies. 
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2 Background on missions 

This section covers the science goals and technical overview of the three missions. 

2.1 WFIRST 

WFIRST is the top-ranked large space mission of the 2010 New Worlds New Horizon (NWNH) Decadal 
Survey. With the addition of a coronagraph, WFIRST would also satisfy the top ranked medium space 
priority of NWNH. The mission is designed to settle essential questions in dark energy, exoplanets, and 
infrared astrophysics. The mission will feature strategic key science programs plus a vigorous program of 
guest observations. A Science Definition Team and a Study Office at the Goddard Space Flight Center and 
Jet Propulsion Laboratory are studying the mission with reports^ in 2013, 2014, and 2015. 

WFIRST is now baselined with an existing 2.4 m telescope NASA acquired from the National Recon¬ 
naissance Office (NRO), a telescope that became available after the completion of NWNH. This configu¬ 
ration is referred to as WFIRST-AFTA (Astrophysics Focused Telescope Asset). The mission consists of 
telescope, spacecraft, a Wide-Field Instrument (WFI), an IFU spectrometer and a coronagraph instrument, 
which includes an IFS spectrometer. The WFI operates with multiple bands covering 0.7 to 2.0 micron 
band with extension to 2.4 microns under study. It has a filter wheel for multiband imaging and a grism for 
wide-held spectroscopy (R=550-800). The pixel scale is 0.11 arcsec, which fully samples the PSF at the H 
band. The IFU has a 3 arcsec FoV and R=75 resolution. The coronagraph operates in the visible 400 - 1000 
nm band. It has a 2.5 arcsec FoV, 10“® effective contrast and 100-200 mas inner working angle. The IFS 
has R=70 resolution. 

WFIRST will measure the expansion history of the Universe as a function of cosmic time and the growth 
rate of the large-scale structure of the Universe as a function of time to test the theory of general relativity 
on cosmological scales and to probe the nature of dark energy. It will employ hve different techniques: su¬ 
pernovae, weak lensing, baryon acoustic oscillations (BAG), redshift space distortions (RSD), and clusters. 
With its 2.4-meter primary mirror, the mission will measure more than double the surface density of galaxies 
detected by the Euclid mission and push these measurements further in the NIR. This higher density will 
enable more detailed maps of the dark matter distribution, measurements not only of two-point statistics 
but also higher order statistics, and multiple tests of systematic dependance of cosmological parameters on 
galaxy properties. The NRC noted in its recent review^: “For each of the cosmological probes described in 
NWNH, WFIRST/AFTA exceeds the goals set out in NWNH.” 

For exoplanets, WFIRST will complete the census of exoplanets begun by Kepler. It will make mi- 
crolensing observations of the galactic bulge using the WFI to discover several thousand planets at the orbit 
of Earth and beyond and will be sensitive to planets as small as Mars. WEIRST’s microlensing survey will 
have the unique capability to detect free-floating planets, thus, enabling astronomers to determine the effi¬ 
ciency of planet formation. The coronagraph will directly image and characterize tens of Uranus to Jupiter 
mass planets around nearby stars and study debris disks. 

WEIRST will also be a worthy successor to the Hubble Space Telescope. With 200 times the field-of- 
view of HST, and the same size primary mirror, it will conduct a rich program of general astrophysics and 
reveal new insights and discoveries on scales ranging from the nearest stars to the most distant galaxies. 

Table 1 gives the WEIRST capabilities for surveys and SN monitoring observations. 

High Latitude Survey (HLS): The nominal HES will be for 2 years with 1.3 years for imaging and 0.6 
years for grism spectroscopy. The coverage will be 2200 deg^ of high Galactic latitude sky within the ESST 
footprint. The imaging will have 2 passes over the survey footprint in each of the 4 imaging filters (J, H, 
El 84 [for shapes] and Y [for photo-z’s]). Data from ESST will be required to provide optical filters for 

^See SDT web site: http://WFlRST.gsfc.nasa.gov/add/ 

http://www.nap.edu/catalog.php7record Jd=l 8712 
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Attributes 

WFIRST Capability 

Imaging survey 

J ~ 27 AB over 2400 sq deg 


J ~ 29 AB over 3 sq deg deep fields 

Multi-filler photometry 

Filters: z, Y, J, H, El84 (Ks), W (wide) 

Slitless wide-field speelroseopy 

0.28 sq deg, R~600 

Slit multi-field speelroseopy 

IFU, R~70 

Number of SN la 

2 X 10^ to z~1.7 

Number galaxies with speetra 

2 X 10^ 

Number galaxies with shapes 

4 X 10« 

Number of galaxies deteeted 

few xlO'^ 

Number of massive elusters 

4 X 10^ 


Table 1: WFIRST capabilities in a nominal ~ 2.5 year dark energy survey 

photo-z determination. Eaeh pass will inelude four 184 see exposures in eaeh filter (with five exposures 
in J band) wifh eaeh exposure offsef diagonally by slighfly more fhan a defeefor ehip gap. This paffern is 
repealed aeross Ihe sky. The speelroseopie survey will have four passes lolal over fhe survey foolprinl wifh 
Iwo “leading” passes and Iwo “Irailing” passes lo enable fhe single grism lo rolale relative lo fhe sky. Eaeh 
pass ineludes Iwo 362 see exposures wifh offsel lo eover ehip gaps. The Iwo “leading” passes (and Iwo 
“Irailing”) are rofaled from eaeh olher by ~ 5°. If is anlieipaled lhal fhe sky eoverage of fhe HES will be 
signifieanlly expanded in fhe extended phase of fhe mission afler fhe baseline firsl 6 years. WEIRST has no 
eonsumables lhal would prohibil a mission of 10 years or longer, and is being designed wifh servieeabilily 
in mind lo enable Ihe possibility of an even longer mission. 

Supernova Survey: A six monlh SN la survey will employ a Ihree-lier slralegy so lo Iraek supernova 
over a wide range of redshifls: 

• Tier 1 for z<0.4: 27.44 deg^ Y=27.1, J= 27.5 

• Tier 2 for z<0.8: 8.96 deg^ J=27.6, H= 28.1 

• Tier 3 for z<1.7: 5.04 deg^ J=29.3, H=29.4 

Tier 3 is eonlained in Tier 2 and Tier 2 is eonlained in Tier 1. Eaeh of Ihese fields will be visited every 5 
days over Ihe 6 monlhs of Ihe SN survey. The imager is used for SN diseovery and Ihe lEU speelromeler is 
used lo determine SN type, measure redshifls, and oblain lighleurves. Eaeh sel of observations will lake a 
lolal of 30 hours of eombined imaging and speelroscopy. The fields are loealed in low dusl regions < 20° 
off an eeliplie pole. A final revisil for eaeh largel for speelroseopy will oeeur after Ihe SN fades for galaxy 
sublraelion. 

Guest Observer Program: A signifieanl amounl of observing time will be awarded lo Ihe eommunily 
Ihrough a peer-seleeled GO program. An example observing program prepared by Ihe SDT has 25% of 
Ihe baseline 6 years of Ihe mission, or 1.5 years, for guesl observations. Signifieanlly more time would be 
awarded in an extended phase after Ihe 6lh year. The GO program is expeeled lo eover broad areas of seienee 
from Ihe solar system lo Galaelie sludies lo galaxies lo eosmology. The April 2013 SDT sludy eonlains a rieh 
sel of ~50 potential GO seienee programs lhal are uniquely enabled by WEIRST. Aslronomieal eommunily 
members provided a wide-ranging sel of one-page deseriplions of differenl GO programs lhal highlighl Ihe 
Iremendous potential of WEIRST lo advanee many of Ihe key seienee questions formulated by Ihe Deeadal 
survey. An imporlanl eonlribulion lhal would likely eome from key projeels in Ihe GO program would be lo 
have deep IR observations of seleeled few-square-degree fields. These ullra deep drilling fields eould reaeh 
limils of J=30 AB. 
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2.2 LSST 


The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-held, ground-based facility designed 
to perform many repeated imaging observations of the entire southern hemisphere in six optical bands (u, 
g, r, i, z, y). The majority of the southern sky will be visited roughly 800 times over the ten-year duration 
of the mission. The resulting database will enable a wide array of diverse scientihc investigations ranging 
from studies of moving objects in the solar system to the structure and evolution of the Universe as a whole 
(Ivezic et al 2014). 

The Observatory will be sited atop Cerro Pachon in Northern Chile, near the Gemini South and SOAR 
telescopes. The telescope incorporates a 3-mirror astigmatic optical design. Incident light is collected by the 
primary, which is an annulus with an outer diameter of 8.4 m, then reflected to a 3.4-m convex secondary, 
onto a 5-m concave tertiary, and finally into three refractive lenses in the camera. The total field of view is 
9.6 deg^ and the effective collecting aperture is 6.6 m in diameter. The design maintains a 0.2 arc-second 
system point spread function (PSF) across the entire spectral range of 320 nm to 1050 nm. The etendue 
(the product of collecting area and field of view) of the system is several times higher than that of any other 
previous facility. 

The telescope mount assembly is a compact, stiff structure with a high fundamental frequency that 
enables fast slew and settle. The camera contains a 3.2 Gigapixel focal plane array, comprised of roughly 
200 4K X 4K CCD sensors with 10 ^m pixels. The sensors are deep depleted, back-illuminated devices 
with a highly segmented architecture that enables the entire array to be read out in 2 s or less. 

Four major science themes have motivated the definition of science requirements for LSST: 

• Taking a census of moving objects in the solar system. 

• Mapping the structure and evolution of the Milky Way. 

• Exploring the transient optical sky. 

• Determining the nature of dark energy and dark matter. 

These four themes stress the system in different ways. For the science of dark energy and dark matter, 
LSST data will enable a variety of complementary analyses, including measurements of cosmic shear power 
spectra, baryon acoustic oscillations, precision photometry of Type la supemovae, measurements of time- 
delays between the multiple images in strong lensing systems, and the statistics of clusters of galaxies. 
Collectively, these will result in substantial improvements in our constraints on the dark energy equation of 
state and the growth of structure in the Universe, among other parameters. 

The main fast-wide-deep survey will require 90% of the observing time and is designed to optimize the 
homogeneity of depth and number of visits. Each visit will comprise two back-to-back 15 second exposures 
and, as often as possible, each field will be observed twice per night, with visits separated by 15-60 minutes. 
Additional survey areas, including a region within 10 degrees of the northern ecliptic, the South Celestial 
Pole, and the Galactic Center will be surveyed with either a subset of the ESST filter complement or with 
fewer observations. The remaining 10% of the observing time will be used on mini-surveys that improve the 
scientific reach of the ESST. Examples of this include the “Deep Drilling Eields” where each field receives 
approximately 40 hour-long sequences of 200 exposures. When all 40 sequences and the main survey visits 
are coadded, this would extend the depth of these fields to r ~ 28. The ESST has identified four distant 
extragalactic survey fields to observe as Deep Drilling Eields: Elias SI, XMM-ESS, Extended Chandra 
Deep Eield-South, and COSMOS. 

Table l^gives the ESST baseline design and survey parameters 

The rapid cadence of the ESST will produce an enormous volume of data, roughly 15 Terabytes per 
night, leading to a total dataset over the ten years of operations of over a hundred Petabytes. Processing 
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Attributes 

LSST Capability 

Einal f-ratio, aperture 

Eield of view, etendue 
Exposure time 

Main survey area 

Pixel count 

Plate scale 

Wavelength coverage 

Single visit depths, design 
Mean number of visits 

Pinal (coadded) depths 

f/L234, 8.4 m 

9.6 deg^, 319 m^deg^ 

15 seconds (two exposures per visit) 
18,000 deg2 

3.2 Gigapix 

50.9 ;um/arcsec (0.2” pix) 

320 - 1050 nm, ugrizy 

23.7, 24.9, 24.4, 24.0, 23.5, 22.6 

56, 80, 184, 184, 160, 160 

25.9, 27.3, 27.2, 26.8, 26.3, 25.4 


Table 2: LSST Capabilities for the fast-wide-deep main survey 


such a large dataset and archiving it in useful form for access by the community has been a major design 
consideration for the project. The data management system is configured in three layers: an infrastructure 
layer consisting of the computing, storage, and networking hardware and system software; a middleware 
layer, which handles distributed processing, data access, user interface, and system operations services; and 
an applications layer, which includes the data pipelines and products and the science data archives. There 
will be both a mountain summit and a base computing facility (located in La Serena, the city nearest the 
telescope site), as well as a central archive facility in the United States. The data will be transported over 
existing high-speed optical fiber links from South America to the U.S. 

The observing strategy for the LSST will be optimized to maximize observing efficiency by minimizing 
slew and other down time and by making appropriate choices of filter bands given the real-time weather 
conditions. The final cadence selection will be undertaken in consultation with the community. A prototype 
simulator has been developed to help evaluate this process, which will be transformed into a sophisticated 
observation scheduler. A prototype fast Monte Carlo optical ray trace code has also been developed to 
simulate real LSST images. This will be further developed to aid in testing science analysis codes. The 
LSST cadence will take into account the various science goals and can also be refined to improve the synergy 
with other datasets. 

LSST anticipates first light in 2020 and the start of the 10-year survey in 2022. It is fully funded for 
construction as of the summer 2014. 

2.3 Euclid 

Euclid is a medium class mission within the European Space Agency’s (ESA) ‘Cosmic Vision’ program^. 
Euclid was formally selected in 2011 and is currently on schedule for launch in 2020. After traveling to the 
E2 Eagrange point and a brief shakeout and calibration period, Euclid will undertake an approximately 6- 
year survey aimed at “mapping the geometry of the dark Universe.” Euclid will pursue four primary science 
objectives (Eaureijs et al 2011): 

• Reach a dark energy EoM > 400 using only weak lensing and galaxy clustering; this roughly corre¬ 
sponds to 1% errors on Wp and Wa of 0.02 and 0.1, respectively. 

• Measure 7, the exponent of the growth factor, with a la precision of < 0.02, sufficient to distinguish 
General Relativity and a wide range of modified-gravity theories. 

^ http://sci. esa. int/cosmic-vision/ 
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• Test the Cold Dark Matter paradigm for hierarchical structure formation, and measure the sum of the 
neutrino masses with a la precision better than 0.03 eV. 

• Constrain the spectral index of the primordial power spectrum, to percent accuracy when com¬ 
bined with Planck, and to probe inflation models by measuring the non-Gaussianity of initial condi¬ 
tions parameterized by /atl to a la precision of ~2. 

Euclid is optimized for two primary cosmological probes: weak gravitational lensing and galaxy clus¬ 
tering (including BAG and RSD). Therefore, Euclid will measure both the expansion history of the Universe 
and the growth of structure. Euclid comprises a 1.2m Korsch 3 mirror anastigmatic telescope and two pri¬ 
mary science instruments. The visible instrument (VIS) consists of 36 4k x 4k CCDs and will be used to 
take images in a single, wide (riz) filter for high precision weak lensing galaxy shape measurements. The 
light entering the telescope will be split via a dichroic to allow for simultaneous observations with the Near 
Infrared Spectrometer and Photometer (NISP), which will take 3 band imaging (Y, J, H) and grism spec¬ 
troscopy in the 1-2 ^um wavelength range. The NISP imaging is aimed at producing high quality photometric 
redshifts when combined with ground-based optical data from a combination of telescopes in the southern 
and norther hemispheres and the spectroscopy is aimed at producing accurate maps of galaxy clustering 
over 2/3 of the age of the Universe. NISP will contain 16 H2RG 2k x 2k NIR detectors procured and 
characterized by NASA. Together, VIS and NISP will survey the darkest (least obscured by dust) 15,000 
square degrees of the extragalactic sky, providing weak lensing shapes for over 1.5 billion galaxies and emis¬ 
sion line spectra for several tens of millions of galaxies, while taking full advantage of the low systematics 
afforded by a space mission. 

While the Euclid hardware and survey design are specifically optimized for dark energy studies, the 
images and catalogs will enable a wide range of ancillary science in cosmology, galaxy evolution, and other 
areas of astronomy and astrophysics. These data will be made publicly available both in Europe via ESA 
and The Euclid Consortium and in the US via the Euclid NASA Science Center at IPAC (ENSCI) after a 
brief proprietary period. Small areas of the Euclid survey (suitable for testing algorithms and pipelines but 
not large enough for cosmology) will be released at 14, 38, 62, and 74 months after survey operations start 
(these are referred to as “quick releases”). The full survey data will be released in three stages: circa 2022, 
2500 square degrees; circa 2025, 7500 square degrees; circa 2028, 15,000 square degrees. 
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3 Enhanced science 


The main dark energy probes established over the last deeade are deseribed next in the eontext of the three 
missions. They are: Type la supernova, galaxy elustering (baryon aeoustie oseillations and redshift spaee 
distortions), weak lensing, strong lensing and galaxy clusters. The biggest advantage from combining infor¬ 
mation from the three missions is in the mitigation of systematic errors, especially via redshift information. 
For the observational and theoretical issues underlying the discussion here we refer the reader to recent 
review articles in the literature (Weinberg et al 2013; Joyce et al 2014). 

LSST, WFIRST and Euclid each have their own systematic errors, most of which are significantly re¬ 
duced via the combination of the survey datasets. The systematic errors affecting each of the surveys indi¬ 
vidually arise from their incomplete wavelength coverage (creating photo-z systematics), their differences 
in imaging resolution and blending (creating shear systematics), or from different biases in galaxy sample 
selections. In most respects the surveys are complementary in the sense that one survey reduces or nearly 
eliminates a systematic in the other. To achieve this synergy requires some level of data sharing; we discuss 
below the cases where catalog level sharing of data is sufficient and others where it is essential to jointly 
process the pixel data from the surveys. 

This section describes each dark energy probe, with a focus on mitigation of systematics, and returns 
to photometric redshifts (photo-z’s). The last three subsections are devoted to stellar science, galaxy and 
quasar science, and the variable universe. 

3.1 Weak lensing 

A cosmological weak lensing analysis in any of these surveys will use tomography, which involves dividing 
the galaxy sample into redshift slices and measuring both the auto-correlations of galaxies within slices and 
cross-correlations of galaxies in different slices. Tomography gives additional information about structure 
growth beyond a strictly 2D shear analysis. However, there are a number of requirements on the data for a 
tomographic weak lensing analysis to be successful. First, we require a good understanding of both additive 
and multiplicative bias in the shear estimates, including how they scale with redshift. Second, we require a 
strong understanding of photometric redshift bias and scatter. This is true in general, but it becomes even 
more important when trying to model theoretical uncertainties such as the intrinsic alignments of galaxy 
shapes with large-scale structure. 

There are a number of ways in which the weak lensing science from each of LSST, WFIRST and Euclid 
would benefit by sharing data with the other surveys and thus improving our understanding of shear estimates 
and/or the photometric redshift estimates. This data sharing may include catalog-level data or possibly image 
data. Although the latter would pose additional complications in terms of the data reprocessing required to 
appropriately incorporate both ground and space images, the gains may be worth the effort where confusion 
or depth or wavelength limitations are severe. 

The main benefit to LSST from either WEIRST or Euclid will be in terms of the calibration of shear 
estimates. The space-based imaging will have much better resolution due to the much smaller PSE size. The 
shapes of galaxies measured from space images will suffer less from model bias and noise bias (at a fixed 
depfh) fhan fhe shapes of fhe same galaxies measured from ground images. LSST does go deeper over ifs 
wide survey fhan fhe space missions, so if will do heller on fhe ouler parls of galaxy images. Eor subsels of 
galaxies measured al comparable depfh, a comparison of fhe shear signal eslimaled by bolh LSST and eilher 
Euclid or WEIRST can help quantify possible syslemalic errors in fhe LSST eslimales. While each survey 
will have ifs own independenl scheme lo correcl for calibralion biases, a cross-comparison of Ihis sorl can be 
used lo validale Ihose schemes. Any addilional correcfion lhal is derived could Ihen be extended fo fhe resl 
of fhe LSST area where Ihere is no such overlap. Nole lhal such a comparison would use a malched sample 
of galaxies in LSST and one of fhe space-based surveys, bul would nol involve comparing per-galaxy shape 



estimates. With different effective resolution, wavelength, and possibly shear estimation methods, there is 
no reason to expect agreement in per galaxy shears (even allowing for variation due to noise). Instead, such 
a comparison would utilize the reconstructed lensing shear fields from the different surveys. WFIRST may 
be more effective than Euclid for this purpose because its narrower passbands mean that it does not suffer 
from the wide-hand chromatic PSF issue mentioned below. 

Another benefit to LSST relates to blended galaxies. The higher resolution space-based images make it 
easier to reliably identify blended galaxies than in ground-based images. The performance of the deblending 
algorithm in LSST image processing pipeline can be significantly improved by providing higher-resolution 
space-based data because of order half of all observed galaxies (though a smaller fraction of the LSST 
gold sample) will be significantly blended with another galaxy. With space-based catalogs, one can at least 
identify which objects are really multiple galaxies or galaxy-star blends and do forced fitting using the 
space-estimated positions. Euclid will be more effective than WFIRST for this purpose because its footprint 
has a greater overlap with that of LSST. 

One benefit to Euclid from combination with LSST relates to chromatic effects. The diffraction-limited 
PSF is wavelength dependent and therefore differs for stars and galaxies due to their different SEDs. Euclid’s 
very wide band imaging means that with its optical imaging alone, it will not be able to correctly estimate 
the appropriate PSF for each galaxy. Multi-band photometry is required to correct for this effect. Without 
correctly accounting for the color-dependent PSF, the systematic will dominate Euclid’s WE error budget, 
so LSST photometry will be very useful in reducing this systematic error. Note that LSST cannot help with 
chromatic effects involving color gradients within galaxies, for which higher resolution imaging data will 
be necessary to estimate corrections; but LSST can provide information about chromatic effects involving 
the average galaxy colors. 

Finally, there are some shear systematics tests that can be done in a joint analysis but not in any one 
survey individually. For example, one route that has been proposed within surveys to mitigate additive 
shear systematics is to cross-correlate shear estimates in different bands or in different exposures rather than 
auto-correlate shear estimates. However, some sources of additive systematics may persist across bands or 
exposures, which would make them quite difficult to identify within the data from the survey alone. Indeed, 
in essentially every survey that has been used for weak lensing analysis to date, the preliminary results 
indicated some unforeseen systematics. Despite care taken in the design and planned operations of new 
surveys, we cannot reliably exclude the possibility of unforeseen systematics. 

In the case of these three surveys, the systematic errors from each survey are expected to be very differ¬ 
ent. Thus, cross-correlating the shear estimates from one survey with those of another should remove nearly 
all shear systematic error contributions, leaving only the true weak lensing signature. If this cross-correlation 
analysis reveals differences from the analysis within individual surveys, it could signal a problem that needs 
to be investigated and mitigated. Furthermore, the surveys have different redshift distributions for the lensed 
galaxies, so a joint cosmological analysis using data from two or more of them will be complementary. 
The joint analysis would also be better able to constrain nuisance parameters like intrinsic alignments and 
photometric redshift errors. 

As described below, the principal way in which Euclid and WFIRST would benefit from sharing infor¬ 
mation with LSST is in the photometric redshift catalog. Without this additional data, Euclid and WFIRST 
will not be able to complete their goals with weak lensing tomography. LSST will also gain from the IR 
band photometry of WFIRST and Euclid in terms of photometric redshift determination. 

The easiest form of data sharing between surveys is if only catalog information is shared. For some 
of the gains mentioned above, this level of sharing is completely sufficient: cross-correlations and shear 
calibration, for example. The deblending improvements would require an iterative reprocessing of the data, 
using a space catalog as a prior and then reprocessing the LSST data with that extra information. A similar 
iterative reprocessing would be necessary for Euclid to properly account for its color-dependent PSF using 
the LSST SED estimates as a prior. For the photometric redshift improvement, one could just use the 
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multiple catalogs, but to do it properly, it is important to have access to the pixel data from both experiments 
at once to properly estimate robust colors for each galaxy. It will be a significant technical challenge to 
develop a framework that will work with the different kinds of image data in a coherent analysis. 

3.2 Large-scale structure 


■ LSST ■ Euclid ■ WFIRST 

■ WFIRST+LSST+Euclid 



Cr(fT8) ^ ^ 


Figure 1! The chart shows how the complementarity of LSST, Euclid and WFIRST contributes to significant improvement in constraints 
on cosmological parameters. As described in the text, the improved constraints on erg come from the mitigation of intrinsic alignment and other 
systematics in weak lensing; the improved constraints on the sum of neutrino masses rriu (in eV) comes from the combination of the weak 
lensing, CMB convergence maps, and galaxy clustering, in particular by reducing the multiplicative bias in shear measurement. Note that the space 
based surveys are assumed to have used ground based photometry to obtain photo-z’s. 


Measurements of large-scale structure probe dark energy via baryonic acoustic oscillation features in 
the galaxy power spectrum. In addition, the full shape and amplitude of the power spectrum measure the 
clustering of matter, if one can also measure the biasing of galaxies relative to the matter. This provides new 
avenues to measure the properties of dark energy and probe gravity. 

The upcoming generation of spectroscopic surveys will make detailed maps of the large-scale structure 
of the universe. DESI, PFS, Euclid and WEIRST will focus on measuring galaxy clustering by obtaining 
spectra of tens of millions of galaxies over redshift ranges from 2 ~ 0.4 — 3.5. The surveys complement 
each other using a variety of spectral lines (Ha, Oil, CII) and galaxy types (EEG and ERGs) as tracers and 
through having distinct, but overlapping, redshift ranges. DESI and Euclid will be wider surveys, while PES 
and WEIRST will go deeper on smaller areas of the sky. 

The combination of spectroscopic and lensing data enables a test of general relativity on cosmological 
scales. The spectroscopic data provides information about galaxy positions and motions determined inertial 
masses moving in the local gravitational potential. The imaging data, by contrast, shows the effect of 
space-time curvature on the trajectories of photons. If general relativity is the correct description of gravity 
then these will agree. Spectroscopic and photometric galaxy cluster surveys from WFIRST, ESST, Euclid 
and DESI, and others, and their cross-correlation with contemporary CMB polarization (CMB lensing) and 
temperature (CMB lensing and kinetic Sunyaev-ZeTdovich) data will provide coincident, complementary 
dynamical and weak lensing tests of gravity in halos two to three orders of magnitude more massive than 
the galaxies within them. Thus we expect powerful tests of gravity by a joint analysis of the three surveys. 
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We study the complementarity of WFIRST, LSST and Euclid by forecasting constraints on matter clus¬ 
tering (we pick CJ8 for convenience, which is sensitive to dark energy and is partially degenerate with the 
parameters, n^, Qm), and the sum of neutrino masses in the Universe. While we concentrate on these 
particular parameters, there are many other possibilities considered by the community. 

Figure [T] shows constraints on as and the sum of neutrino masses from the three surveys individ¬ 
ually and in combination. We include a CMB prior from Planck, and marginalize over cosmological (and 
galaxy bias) parameters when appropriate. For the space surveys we assume that photo-z’s will be obtained 
with the aid of ground based data. WFIRST’s densely sampled spectroscopic galaxies will provide a de¬ 
tailed characterization of filamentary structure that is complementary to the lensing measurements provided 
by imaging surveys. This allows us to understand one of the most elusive systematics in weak lensing sur¬ 
veys - the intrinsic alignment (lA) of galaxies. One can test and validate lA models and marginalize over 
the free parameters in specific models when we can have a 3D map of the universe via galaxy spectra. In 
FSST, the lensing sources spans a large redshift range between 0 and 3 and a large luminosity range. With 
WFIRST and Euclid, we will have spectroscopic sources spanning a large part of this redshift range, which 
enables a very good construction of filament map of the surveyed volume. WFIRST offers a much higher 
density of sources while FSST and Euclid will cover a much large area of sky. The joint analysis leads to 
significant improvement in the constraint on fjg, as shown on the left in Figure [T] 

We also consider the constraints on massive neutrinos which suppress the growth of structure below the 
free-streaming scale. The combination of WFIRST and Euclid’s spectroscopic sample, along with CMB 
lensing convergence, can be used to calibrate the shear multiplicative bias in FSST sources. We follow 
Das, Evrard & Spergel (2013) in methodology and calculate the improvement in constraints first in the 
multiplicative bias in shear, and propagate this to the constraints on the sum of neutrino masses that relies 
on the large sky coverage of FSST. We assume that a Stage III level CMB convergence field will be available. 
We plot the improvement in Tiniu on the right in Fig[T] 

3.3 Galaxy clusters 

FSST, Euclid and WFIRST will use galaxy clusters (primarily via gravitational lensing mass measurements) 
as part of their cosmological measurements. However, both for cosmology and the study of galaxy evolution 
in clusters, a joint analysis of the data can greatly improve the reach of each mission. Fundamentally, this is 
because each of the missions observes different aspects of the emission from clusters. FSST provides optical 
colors, useful for cluster identification via photometric redshifts, as well as (for group-scale structures) time 
delays for lensed objects. Euclid provides shallow NIR photometry and somewhat deeper NIR spectroscopy, 
as well as higher resolution optical imaging. WFIRST will provide much better resolution images and deep 
NIR photometry and spectroscopy over a smaller area than Euclid. Because of the complementarity of 
these measurements, combining the data has implications for cluster finding, redshift determination, mass 
measurement and calibration, and the study of galaxy evolution in clusters. 

Galaxy cluster finding: Over the past decade, techniques have been developed that select galaxy groups 
and clusters very effectively by isolating overdensities simultaneously in projected position and color space 
(for example redMaPPer). These techniques have been optimized for optical wavelengths. FSST will be ex¬ 
tremely efficient at finding groups and clusters of galaxies and estimating their redshifts (although eROSITA 
will also be an efficient cluster finder, the redshift information available from the FSST imaging will allow 
selection of samples at fixed redshift for study). Although redMaPPer is effective to z ~ 1, extending its 
redshift range will require infrared imaging, which would be provided by Euclid or WFIRST. Given that the 
number density of detected (as opposed to resolved) sources in FSST and WFIRST are going to be similar, 
and given the importance of consistent photometry measurements across bands, simultaneous measurements 
of the pixel data would yield more accurate selection of clusters. 
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Redshift Determination: As discussed below, the photo-z’s of galaxies will be greatly improved by com¬ 
bining the optical and infrared imaging of the three missions. Thus, joint analysis of the LSST-i-WFIRST/Euclid 
data gives more accurate cluster sample selection, and allows accurate determination of the background 
galaxy redshift distribution, which is the largest single source of systematic uncertainty in determination of 
the mass function from weak lensing mass determinations. 

Strong and Weak Lensing: Measuring the weak lensing shear requires high resolution imaging. WFIRST 
will be able to resolve many more galaxies than LSST, and as discussed above will aid the shear calibration 
of LSST. In addition, the conversion of shear to mass depends on redshift information. In this way, the 
combination of LSST, Euclid, and WFIRST data will provide much more accurate mass measurements than 
either will alone. 

Cluster arc tomography is greatly improved by the joint pixel-level analysis of LSST, Euclid and WEIRST 
data. Eor instance, WEIRST will have the angular resolution to discover ~ 2000 strong lensing clusters in 
its survey footprint. However, most of the arcs will be too faint (and too low surface brightness) to have 
spectroscopic redshifts. By contrast, ESST will reach surface brightness limits (> 28.7 Mag per square 
arcsecond) sufficient to detect the arcs and obtain photometric redshift estimates, but will not have the res¬ 
olution to separate the arcs cleanly from the cluster galaxies in many cases. Thus, the measurement of the 
tomographic signal from multiple arc systems, which is a strong test of cosmology, will be much stronger 
from the joint analysis than from Euclid, WEIRST or ESST alone. Once again, the added value of the pixel 
data is great, because cluster cores are very dense environments and obtaining valid photometric redshifts 
will require use of the space-based imaging to disentangle the light contribution form the arcs and cluster 
galaxies. 

Systematics: One of the sources of systematic uncertainty in weak lensing shear measurements from 
galaxy clusters comes from the blending of sources both with other sources and with cluster galaxies. The 
source-source blending is in many ways similar to that encountered in large-scale weak lensing measure¬ 
ments, but in clusters, there is an additional source of bias from blending with cluster galaxies. This is 
both radially dependent and (because cluster concentration evolves with redshift) redshift dependent, and 
directly affects the normalization of the cluster mass function. Using the WEIRST data to directly measure 
the deblending bias in the clusters in common between the two surveys would make it possible to correct 
for the bias. 

3.4 Supernovae 

Each individual Type la supernova (SN la) is potentially a powerful probe of relative distance in the Uni¬ 
verse. Thousands of them provide detailed information about the expansion history of the Universe, in 
particular allowing us to probe the accelerated expansion and the nature of dark energy. Hundreds of thou¬ 
sands of them will allow us to probe the nature of dark energy in different regions and environments. A set 
spanning the past 10 billion years of cosmic history will allow for a rich and powerful exploration of the 
nature of the expansion and more recent acceleration of the Universe and thus reveal key insights into the 
key constituents. 

Ground-based SN la searches have provided almost one thousand well-observed SNe la in 20 years of 
work. Increasingly powerful ground facilities will rapidly increase this number to thousands per year (e.g., 
Pan-STARRS, PTE and DES) and then hundreds of thousands per year in the era of ESST. Eor a smaller 
subset (hundreds) of these supernovae, even calibrated spectrophotometric time-series have been obtained 
(at the nearby SN Eactory). However, the rest-frame emission of SNe la is line-blanketed at UV wavelengths. 
Thus, by z ~ 1 their light has already redshifted out of the optical atmospheric window easily accessible 
from the ground. In addition, because the relative distances across redshifts are estimated based on relative 
flux differences across wavelength, accurate cross-wavelength calibration is critical to use supernovae to 
determine the properties of dark energy. 
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The stability of WFIRST in space yields significant benefits in both wavelength coverage and flux cal- 
ibrafion for fhousands of SNe la/year from 0.1 < 2 : < 1.7. The accessibilify of fhe ground allows for 
fhe large aperfure and high-efendue of fhe LSST sysfem fhaf can observe 10,000s of SNe la per year from 
0.03 < z < 1. 

Measuring disfances from 2 ; ~ 0.03 fo 2 ; ~ 1.7 will conned fhe era when dark energy firsf made ifs 
presence known, fhrough fhe fransifion fo an accelerafing expansion, and up fo fhe presenf day. These are fhe 
key redshiffs over which we will discover hinfs fo fhe underlying nafure of dark energy, and fhe combinafion 
of fhe high-redshifl reach and calibrafion of WFIRST and fhe huge numbers and volume coverage of LSST 
will yield a sfrong and powerful pairing for SN Cosmology. 

The key question for Ibis work will be fhe confrol of sysfemafic uncerfainlies as we explore a large 
span of time and changing populafions of hosf-galaxy environmenfs in which fhe supernovae are born. 
The combinafion of WFIRST, LSST, and ofher ground-based follow-up insfrumenfs will make possible 
unprecedenfed rich characferizafion of each supernova so fhaf low- and high-redshifl supernovae can be 
malched al a more delailed level fhan wifh fhe broad “Type la” definilion. For example, a lime series of 
delailed speclropholomefry could be oblained for fhe same resl-wavelengfhs for every supernova. 

Modern SN dislance analyses are already joinl analyses. The currenl Hubble Diagrams are conslrucled 
from supernovae from differenl surveys (e.g., Befoule el al. 2014 combined fhe Supernova Legacy Survey 
(SNLS), SDSS Supernova Survey, and fhe currenl heterogeneous sample of nearby 2 : < 0.05 SNela). Buf 
a frue joinl analysis should begin al fhe pixel level and lake advanlage of all fhe internal delails aboul 
calibrafion and pholomelry wilhin each survey. Specifically, when combining SNela from differenl dafa 
sources, a joinl analysis is beneficial because of: 1. Consislenl pixel-level analysis, in particular pholomelric 
exlraclion; and, 2. Deep cross-checks of pholomelric calibration nol possible wifh jusl exlracled SN lighf 
curves. 

Supernova surveys have dislincl science requiremenls bolh on supernova discovery and early classifi- 
calion, and on fhe qualify of dislance modulus and redshifl delerminalions of discovered supernovae. Each 
science requiremenl propagates info hardware and survey requiremenls. Technical requiremenls for super¬ 
nova discovery (e.g., large field of view, detection of new poinl sources) can have a more cosl-effeclive 
ground-based solution, whereas requiremenls for pholomelry and speclroscopy (e.g., high signal-lo-noise 
SN speclral fealures, infrared wavelenglh coverage) are only fulfilled wilh a space observatory. A case in 
poinl is Ihe baseline WFIRST plan where an imaging survey is conducted to identify supernovae for largeled 
speclroscopic observations. The large field of view of LSST and Ihe availabilily of SN photons deleclable 
by LSST oul to 2 : = 1.2 should make Ihe LSST Ihe belter inslrumenl for SN discovery, especially al low 
redshifls where Ihe angular density of potential hosl galaxies is low. Assuming Ihe Iriggering logistics work. 
Ibis mighl free up WFIRST to do whal if is belter al: oblaining supernova speclropholomelric time series. 

In addition, Ihere is a benefil of observing Ihe same SN wilh multiple overlapping surveys. Such obser¬ 
vations provide: 

• Direcl cross-inslrumenlal calibration on Ihe source of interest (If is imporlanl to note, however, lhal 
Ibis calibration is nol Ihe same as Ihe cross-wavelenglh calibration lhal is required when comparing 
supernovae al very differenl wavelenglh.) 

• Expanding observations of Ihe wavelenglh/lemporal range of Ihe supernova. 

• A quanlilalive assessmenl of systematic errors. 

Einally, supernova cosmology analyses require a suite of external additional resources beyond Ihe direcl 
observations of high-redshifl supernovae. Thus, planning for a joinl ESST-WEIRST SN cosmology analysis 
will yield signilicanl gains in combining resources to oblain external dala and Iheory needs. Some key 
examples include: 
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• Low-z SN set, including possible spectrophotometric time-series follow up. 

• Fundamental flux calibrators (which can provide the cross-wavelength calibration that is required 
when comparing supernovae at very different wavelength). 

• High-resolution spectroscopy of host galaxies. 

• Improved SN la theory and empirical modeling. 

3.5 Strong tensing 

The LSST, Euclid and WFIRST datasets will be particularly complementary in the field of strong gravita¬ 
tional lensing. Galaxy-scale strong lenses can be used as probes of dark energy (via time delay distances 
or multiple source plane cosmography). The time delays between multiple images in a strong gravitational 
lens system depend, as a result of the lens geometry, on the underlying cosmology - primarily the Hubble 
constant but also, in large samples where internal degeneracy breaking is possible, the dark energy parame¬ 
ters. (Refsdal 1964, Tewes et al 2013; Suyu et al 2014). Lens systems with sources at multiple redshifts may 
also provide interesting and competitive constraints on cosmological parameters, via the ratio of distance 
ratios in each system. 

Strong lenses also provide unique information about dark matter, via the central density profiles of 
halos on a range of scales, and the mass function of sub-galactic mass structures that cause measurable 
perturbations of well-resolved arcs and Einstein Rings. Both this science case and strong lensing cosmology 
will benefit greatly from the large samples of lenses detectable in these wide field imaging surveys; we sketch 
these cases out below in Section 4, and identify how high fidelity space-based imaging and spectroscopy 
will enable new measurements. 

3.6 Photometric redshifts 

Photometry provides an efficient way to estimate the physical properties of galaxies (i.e. their redshifts, 
spectral types, and luminosities) from a small set of observable parameters (e.g. magnitudes, colors, sizes, 
and clustering). LSST will depend on these photometric redshifts for all of its major probes of dark energy 
as it is infeasible to obtain spectroscopic redshifts for the bulk of objects in the LSST sample with any 
reasonable amount of telescope time. For the WFIRST and Euclid missions, three of the key applications 
will rely on photometric redshifts: measures of the mass power spectrum from analyses of the weak lensing 
of faint galaxies, breaking of the degeneracies between possible redshift solutions for galaxies exhibiting 
only a single emission line in grism spectroscopy, and the identification of high redshift Type la supemovae 
based on the photometric redshifts of their host galaxies. In addition, photometric redshifts will play a key 
role in much of the LSST, WFIRST, and Euclid extragalactic science where we wish to constrain how the 
demographics of galaxies change over time. 

The high efficiency of photometric redshifts comes at a cost. The primary features that enable photo¬ 
metric redshift estimation are the transitions of the Balmer (3650 A) and Lyman (912 A) breaks through a 
series of photometric filters. Uncertainties in isolating the position of the breaks due to the low spectral res¬ 
olution of the photometry give rise to a scatter between photometric redshift estimates and the true redshifts 
of galaxies. Confusion between breaks (due to sparse spectral coverage and incomplete knowledge of the 
underlying spectral energy distributions) leads to catastrophic photo-z errors, as seen in Figure]^ 

Existing studies have demonstrated that while photometric redshift accuracies of cj^ ~ 0.007(1 4- z) 
are possible for bright objects with many filters, uncertainties of Uz ~ 0.05(1 -|- z) are more typical when 
restricted to 5-6 bands of deep imaging. Catastrophic photometric redshift errors (where the difference 
between the photometric and spectroscopic redshifts, Az, exceeds the 3 — a statistical error) typically 
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occur in well over 1% of cases. As the photo-z scatter and catastrophic failure rate increase, information 
is degraded and dark energy constraints will weaken. Furthermore, if there are nonnegligible (>~ 0.2%) 
systematic offsets in photo-z’s or if the Az distribution is mischaracterized, dark energy inference will be 
biased at levels comparable to or greater than expected random errors. As a result, careful calibration and 
validation of the photometric redshifts will be necessary for the LSST, WFIRST, and Euclid surveys. 

3.6.1 The impact of Euclid and WFIRST near-infrared data on LSST photometric redshifts 

The LSST filter system covers the u, g, r, i, z and y passbands, providing substantial leverage for redshift 
estimation from z = 0 to z > 6. For the LSST “gold” sample of galaxies, i < 25.3, Figure shows how 
the HLS for WFIRST, with four bands comparable in depth to the 10 year LSST survey, (i.e. 5 it extended 
source depths of Y^b= 25.6, Jy4s=25.7, Hab= 25.7, and F184 ab= 25.2), will significantly improve on the 
photometric redshift performance from LSST alone. For example, at redshifts z > 1.5, where the Balmer 
break transitions out of the LSST y band and into the WFIRST and Euclid infrared bands, the inclusion 
of the WEIRST data results in a reduction in by a factor of more than two (1.5 < z < 3), and a 
reduction in the fraction of catastrophic outliers to <2% across the full redshift range. Euclid’s three-band 
NIR photometry, while shallower, will have a much greater overlap with ESST and will also provide a 
quantitative improvement in ESST photo-z’s. Combining the ESST, WEIRST and Euclid photometric data 
effectively will depend, however, on the details of the respective filter systems, their signal-to-noise, our 
ability to extract unbiased photometric measurements from extended sources (e.g. the deblending of sources 
using the higher spatial resolution of the WEIRST data), and the accuracy of the photometric calibration of 
the data both across the sky and between the near-infrared and optical passbands. 



Eigure 2: A comparison of the relative photometric redshift performance of the LSST optical filters (left panel) with a combination 
of LSST and WFIRST filters (right panel). The simulated data assumes a 10-year LSST survey and a “gold sample” with i < 25.3. 
The addition of high signal-to-noise infrared data from WFIRST reduces the scatter in the photometric redshifts by roughly a 
factor of two (at redshifts 2 > 1.5) and the number of catastrophic outliers by a factor of three. These simulations do not account 
for deblending errors or photometric calibration uncertainties, and assume complete knowledge of the underlying spectral energy 
distributions of galaxies as an ensemble. 
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Figure 3: Predictions of the fraction of LSST weak lensing sample objects that would yield a secure (multiple-confirmed-feature) 
spectroscopic redshift, based either on 1440-second exposure time with WFIRST (colored regions) or 10 nights’ open-shutter- 
time spectroscopy with the Subaru/PFS spectrograph (black curve) WFIRST IFU spectroscopy would provide training redshifts for 
objects at higher 2 ; than are easily accessible from the ground, particularly if read noise per pixel is small (the colored regions 
indicate a range of feasible scenarios). Longer exposure times (e.g., in supernova fields or by optimized dithering strategies) could 
enhance the success rate further. 

3.6.2 Mitigating systematics with WFIRST and Euclid spectroscopy 

The optimization of photometric redshift algorithms and the calibration of photometric redshift uncertainties 
both require spectroscopic samples of galaxies. If simple algorithms are used, more than 100 spectroscopic 
survey regions (of ~0.25 deg^) with at least 300-400 spectroscopic redshifts per region may be required 
to optimize a photometric redshift algorithm (whether by refining templates and photometric zero points or 
as input for machine learning algorithms) to ensure that their accuracy is not limited by sample variance 
in the spectroscopic training set (Cunha et al. 2012); with techniques that take this variance into account, 
15-30 fields may be sufficient (Newman et al. 2014). An ideal training set would span the full range 
of properties (including redshifts) of the galaxies to which photometric redshifts will be applied. To the 
degree to which we do not meet this goal, we can expect that photometric redshift errors will be degraded, 
weakening constraints on dark energy as well as other extragalactic science. 

Current spectroscopic samples fare well; surveys to i? = 24.1 or i = 22.5 (more than two magnitudes 
shallower than the LSST “gold sample”) with 8-lOm telescopes have obtained >99% secure redshifts for 
21-60% of all targeted galaxies, and >95% secure redshifts for 42-75% of the galaxies (incorrect redshift 
rates above 1% would lead to photo-z systematics that exceed LSST requirements). WFIRST spectroscopy 
can address these limitations in a number of ways. Depending on the final configuration and difhering 
sfrafegy for WFIRST, an objecf in fhe LSST weak lensing sample will fall on fhe IFU field-of-view af leasf 
10% (and up fo ~ 100%) of fhe time. Mosf objecfs in fhis sample would yield a successful redshiff wifh a 
~ 1440 sec exposure (see Figurej^. For an IFU wifh a 3“ x 3“ field-of-view, al leasf 10,000 specfra down 
fo fhe LSST weak lensing deplh, corresponding lo roughly 20,000 down lo fhe WFIRST limils, would be 
measured concurrently with the WFIRST HLS (with minimal impact from sample/cosmic variance). This 
spectroscopy would have very different incompleteness from ground-based samples, allowing a broader 
range of galaxies to have well-trained photometric redshifts, with accuracy limited by the imaging depth 
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rather than our knowledge of galaxy SEDs. 

If our speetroscopic samples do not have a success rate approaching 100%, accurate characterization 
of photo-z errors (required for dark energy inference) is likely to instead be based on the cross-correlation 
between spectroscopic and photometric samples (cf. Newman 2008). To meet LSST calibration require¬ 
ments, such an analysis will require secure (< 1% incorrect) spectroscopic redshifts for ~ 100, 000 galaxies 
spanning the full redshift range of the photometric samples of interest and an area of at least a few hundred 
square degrees. Euclid and WEIRST grism spectroscopy will provide spectra for many millions of galaxies 
over wide areas (e.g., 2 x 10^ galaxies in the WEIRST HES). However, this spectroscopy will provide the 
highly secure, multiple-feature redshifts required for training and calibrating photo-z’s in only very limited 
redshift ranges. As a result, Euclid and WEIRST grism spectroscopy may contribute to photometric redshift 
calibration in combination with other datasets spanning the remaining redshift range (e.g., from DESI), but 
cannot solve this problem on their own. 

3.7 Stellar science 

New insights on stellar populations in the Milky Way and nearby galaxies depend critically on gains in 
photometric sensitivity and resolution. Thus far, the state of the art imaging surveys of our Galaxy have 
involved relatively modest tools with respect to today’s capabilities. Eor example, the SDSS survey is both 
shallow and operates at low resolution, and IR surveys such as 2MASS and WISE are based on even smaller 
telescopes operating at even lower resolution. 

The combination of ESST, Euclid, WEIRST and Gaia will transform our view of the Milky Way and 
nearby galaxies. ESST will map 18,000 square degrees of the night sky every few days, across 6 visible light 
filters. The imaging depth of ESST will be over 100 times (5 magnitudes) deeper than SDSS, and at higher 
resolution. WEIRST will overlap some of the same footprint, but will also extend the spectral range of Milky 
Way surveys to near-infrared wavelengths at even higher spatial resolution (0.1 arcsecond pixels). Euclid 
provides much higher spatial resolution than ESST in the visible (aids deblending sources and crowded field 
phofomefry) and Gaia addresses fhe single biggesf confribufion fo fhe error budge! of mosf sfellar population 
sfudies aimed af characferizing physical properfies; knowledge of fhe disfance of sources. 

Much of fhe sfellar populafion science case for ESST (e.g., from fhe Science Book) overlaps fhaf from 
WEIRST (e.g., fhe communify science whife papers in fhe 2013 WEIRST SDT reporf). These include 
sfudies of fhe inifial mass funcfion, sfar formation hisfories for Milky Way componenfs, near field cosmol¬ 
ogy fhrough deep imaging and proper mofions of dwarf safellifes, discovery and characferizafion of halo 
subsfruclure in nearby galaxies, and more. Given fhe wavelengfh complemenfarify alone, mosf sfellar pop- 
ulafion invesfigafions will be aided by fhe panchromafic baseline of ESST (and Euclid) -i- WEIRST (e.g., 
inferprefing SEDs info fundamenfal sfellar properties, increasing fhe baseline for sfar formafion hisfory fils 
from multi-band color-magnilude diagrams, confirming memberships of very blue or red slars in eilher dala 
sel, elc.). However, fhe overall qualify of all Ihree dafa sels for general invesfigafions can be greally enhanced 
(and, in a uniform way) by joinlly re-processing fhe pixel dafa sefs for ESST, Euclid, and WEIRST. 

To highlighl fhe lechnical requiremenfs for Ihis join! processing, we poinl lo Ihree ESST and WEIRST 
’’killer apps” in sfellar population science, 

• Eslablishing and characferizing fhe complefe speclrum of Milky Way safellifes 

• Testing halo formation models by measuring fhe shapes, subsfruclure conlenl, masses, and ages of a 
large sel of sfellar halos wilhin 50 Mpc 

• Measuring fhe sfar formafion hisfory, Galaclic mass budgel, and spalially and chemically-dependenl 
IME 
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In each of these cases, LSST and Euclid have the advantage with respect to field of view (and cadence). 
WFIRST has the advantage of higher photometric sensitivity to red giant branch and low mass star tracers, 
and much higher resolution for star-galaxy separation in the IR. Joint processing of the three data sets offers 
the following rewards, 

1. By using the WFIRST astrometry to deblend sources and feed new positions of stars to FSST (and 
to a lesser extent, Euclid), the FSST and Euclid photometry can be made to go deeper. In the actual 
processing, the pixels from FSST, Euclid and WFIRST would be analyzed simultaneously to detect 
and perform photometry on all detections in any of the more than dozen FSST-i-Euclid-i-WFIRST 
bandpasses. 

2. WFIRST photometry alone will be critically insensitive to hot stars (e.g., white dwarfs in the Milky 
Way, horizontal branch stars in the Focal Group). In most halo and sparse populations, these objects 
will be easily detected by FSST, and to a shallower depth, Euclid. Similar to 1.), the FSST photometry 
can be used to obtain new measurements in the WFIRST data set and to get an infrared color for these 
sources. 

3. As a result of the above, many more sources in the overlapping footprint of the three telescopes 
will have high-precision panchromatic photometry. This increased wavelength baseline will yield a 
much more accurate characterization of the underlying stellar populations (e.g., the sensitivity to age, 
reddening, and distance variations increases on this baseline). Note, this is different from matching 
FSST, Euclid, and WFIRST catalogs since the goal here is to detect sources that otherwise would not 
have been measured in one of the three data sets. The biggest gain in our interpretation of the stellar 
populations to derive physical properties will come from accurate distances to the stellar populations 
from Gaia. The different photometric sensitivities is not a concern here, since there are bright tracers 
in most populations that Gaia can target. 

4. The WFIRST morphological criteria of faint sources provides a crucial assessment of FSST’s star- 
galaxy separation. 

Many of these joint analyses can be done at the catalog level while others likely require working with the 
pixels from each of the surveys. 

3.8 Galaxy and quasar science 

The WFIRST photometry will go to a similar depth (on an AB scale) as the full coadded FSST images. 
Euclid NIR data, while shallower, will overlap a much larger area of the FSST survey. Thus the combination 
of FSST, Euclid and WFIRST will give high S/N 9-band photometry stretching from 3000 A to 2 microns 
(a factor of six in wavelength, twice that of either WFIRST/Euclid or FSST alone) for the Gold sample 
(i < 25.3) of hundreds of millions of galaxies. This broad wavelength coverage probes the SEDs of galaxies 
beyond 4000A to redshifts 3 and beyond, allowing detailed determination of star formation rates and stellar 
masses during the cosmic epoch when galaxies assembled most of their stars. (Unobscured) AGN will be 
identified both from their SEDs and from their variability in the repeat FSST imaging. 

The Hubble Space Telescope has been able to obtain comparably deep photometry over a similar wave¬ 
length range over tiny areas of sky; for example, the Cosmic Assembly Near-infrared Deep Extragalactic 
Fegacy Survey (CANDEFS) covers only 0.2 deg^, while WFIRST will cover an area ten thousand times 
larger. With this much larger sky coverage, the combination of FSST and WFIRST can explore the evolu¬ 
tion of galaxy properties over this broad range of cosmic history. 

The weak lensing signature from individual galaxies is too small to be measured, but stacking analyses 
can be done in fine bins of stellar mass and redshift, exploring the relationship between dark matter halo 
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mass, stellar populations, galaxy morphology, AGN activity and environments during the peak of galaxy 
assembly. The exquisite resolution of the WFIRST images will also allow the role of mergers in galaxy 
evolution to be quantified, with galaxy pairs separated by as little as 3 kpc discernible over essentially the 
entire redshift range. 

The multi-band photometry will allow both quasars and galaxies at much higher redshifts to be identified 
as well. Currenf sfudies wifh HST have allowed idenlificafion of a handful of galaxy candidates fo redshifl 8. 
WFIRST -I- LSST will expand fhese sfudies over 10,000 limes Ihe area, looking for objecls fhaf are dropoufs 
in fhe u fhrough y fillers. The wide area is parlicularly imporlanl for discovering Ihe rare mosl luminous 
objecls al fhese redshifls, lhal can be followed up in delail wifh JWST and Ihe nexl generation of ~30-meter 
telescopes. The sludy of such objecls probes Ihe exil from Ihe Dark Ages, as Ihe inlergalaclic medium 
became reionized. Underslanding Ihe formation of Ihe firsl galaxies and Ihe reionizalion process in delail 
requires pulling Ihese objecls in Iheir large-scale slruclure conlexl, which is only possible wilh samples 
selected over large solid angles. In particular, Ihe expanding bubbles by which reionizalion is Ihoughl lo 
progress are Ihoughl lo subtend tens of arcminules on Ihe sky, requiring wide-field samples lo explore Ihem. 

3.9 The variable universe 

Time domain aslronomy encompasses an exlremely varied suite of aslrophysical topics, from Ihe repeated 
pulsations of regular variable slars lhal serve as Iracers of Galactic slruclure and cosmological dislances (e.g. 
RR Lyrae and Cepheids), to Iransienl evenls lhal by Iheir nalure occur eilher unprediclably or only once (e.g. 
stellar flares, supernovae, tidal disruption evenls). As suggested by Ihese examples, Iransienls and variables 
occur bolh locally and al Ihe mosl dislanl reaches of our universe, and so comprise a population wilh bolh 
a very high dynamic range in luminosity and a wide span of colors. Time domain aslronomy is a rapidly 
growing field: Ihe rich public dalasel provided by LSST will enable a large pari of Ihe aslronomy community 
to analyze observations of Iransienl and/or variable objecls. The potential synergies belween LSST, Euclid 
and WFIRST are as compelling and varied as time domain aslronomy ilself, bul here we highlighl some 
compelling examples. 

The combined WFIRST, Euclid and ESST dalasels would greally facililale Ihe characterization and 
classification of Iransienl and variable objecls. One of Ihe challenges imminenl in Ihe coming ESST era 
is Ihe identification of samples of variable and Iransienl evenls to be largeled for additional observations. 
Erom Ihe Iransienl side, Ihe anticipated high volume of Ihe ESST alerl slream requires rapid Iriage of large 
numbers of evenls to identify Ihe mosl imporlanl Iransienls: Ihose evenls which, if nol pursued wilh followup 
observations immediately, will be opporlunilies losl forever. Eurlhermore, Ihe high volume of evenls also 
necessilales Ihe prioritization of Ihe mosl compelling sources to followup, which in lurn requires rapid 
characterization of an evenl (where characterization here is distinguished from aclual classification into a 
known category of objecl). The WEIRST and Euclid dalasels will assisl wilh Ihese challenges in a number 
of ways. Eirsl, Ihe deep infrared observations will greally assisl in Ihe identification of Ihe “foreground fog” 
of very red Galactic sources: low mass slars and brown dwarfs. While ESST’s deplh will identify many of 
Ihese objecls in quiescence, some percenlage will escape even ESST on accounl of Iheir very red colors and 
low luminosities. Therefore, flares on Ihese objecls may still create whal appear to be cosmological Iransienl 
evenls in spite of Iheir local origin. The deep, panchromatic characterization of Ihe sialic sky will enable 
Ihe association of Iransienl evenls wilh co-localed sources, whelher Ihose sources are quiescenl detections 
of Ihe same source, or possible hosl galaxies of Ihe evenl whose properties may give clues for a deeper 
underslanding of Ihose Iransienls. 

As pointed oul in Ihe previous subsection on stellar science, Ihe slrenglh of Ihe combined dalasels goes 
far beyond malching Ihe WEIRST, Euclid and ESST calalogues. Ralher if increases Ihe deplh of Ihe pho- 
lomelry by making il possible to delecl sources lhal would olherwise nol have been detected. The ability 
to identify (or rule oul) infrared counlerparls may also be helpful in confirming some of Ihe mosl elusive 
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transients. For example, potential eleetromagnetie eounterparts to gravitational wave sourees are expeeted 
by eurrent models to emit their peak emission in near IR wavelengths. The assoeiation of any gravitational 
wave event detention with a known galaxy greatly reduees the seareh volume for EM eounterparts to grav¬ 
itational wave events, and the elimination of faint red baekground galaxies as putative assoeiated transient 
kilonova emission will help eonfirm potential multi-messenger events. 
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4 Coordinated observations 


A significant fraction of WFIRST time will be devoted to a Guest Observer program. In this section of 
the report, we consider a number of potential GO programs that will likely follow-up on or extend LSST 
observations. While the GO programs will likely be selected by a future Time Allocation Committee (TAG), 
we will benefit by preparing for likely observing scenarios. Here, we consider three examples of joint 
LSSTAVFIRST studies: WFIRST follow-up of strong lenses and lensed AGNs found in LSST, joint obser¬ 
vations of deep drilling fields, and joint WFIRST/LSST searches for the afterglows of gravitational wave 
bursts. 

4.1 Measuring dark energy with strong lens time delays 

As described in the previous section, time delays between multiply imaged sources are a powerful cosmo¬ 
logical probe. LSST will enable time delay distance measurements towards hundreds of lensed AGN and 
supernovae (Oguri & Marshall 2010, LSST Science Book). The detection of candidate systems, which have 
image separations of just 1-3 arcseconds, and confirmation of them by lens modeling, depends critically 
on the quality and depth of the imaging data, suggesting that there is much to be gained by searching a 
combined LSST/EuclidAVFIRST imaging dataset. Measuring cosmological parameters with this sample 
will require high signal to noise ratio and high resolution imaging and spectroscopy of any Einstein Rings 
discovered in order to constrain the lens mass distributions. By integrating longer on the most interesting 
systems, WEIRST can achieve background source densities of over 200 sources/arcmin^, which will enable 
detailed weak lensing mass maps of these systems. 

Joint analysis of all ESST and wide field infrared survey data will improve the completeness and purity 
of the lens searches in the following ways: 1) a combined catalog of ESST and space-detected objects will 
produce a more accurate initial list of target lenses, and 2) joint analysis of all the available cutout images 
will enable better lens modeling, leading to a much cleaner sample of candidates to be confirmed in the 
followup observations. 

Being able to model large proportions of the lens sample precisely with high resolution space-based 
imaging will enable the surveys themselves to constrain the structure and evolution of massive galaxies 
across a homogeneously-selected sample, providing important prior information needed by other cosmology 
probes. 

The joint optical and infrared photometry provided by ESST, Euclid and WEIRST will constrain the 
mass environments of the strong lenses (including the line of sight environment) by providing accurate 
photometric redshifts and stellar masses for all the relevant galaxies in the lightcone (e.g., Greene at al 
2013, Collett et al 2013). 

4.2 Probing dark matter with grism spectroscopy of lensed AGN 

Measuring the mass function of dark matter subhalos down to small masses (e.g., 10® solar masses) is a 
unique test of the CDM paradigm, and of the nature of dark matter in general. The CDM mass function 
is expected to be a power law rising as ® down to Earth-like masses, while a generic feature of self¬ 
interacting dark matter and warm dark matter models is to introduce a lower mass cutoff. Current limits set 
that cutoff at somewhere below 10®M©. We know from observations of the environs of the Milky Way that 
there is a shortage of luminous satellites at those masses, but this may be entirely due to the physics of star 
formation at the formation epochs of these satellites. Dark Milky Way satellites at lower masses may be 
detectable through their dynamical effects on tidal streams, but statistical studies of the one galaxy we live 
in will always carry high uncertainty. 
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A powerful way to detect subhalos independent of their stellar content and in external galaxies is via 
the study of large samples of strong gravitational lenses, in particular the so-called flux ratio anomalies (see 
Treu 2010, and references therein). Small masses located in projection near the four images of quadruply- 
lensed quasars cause a strong distortion of the magnification and therefore alter the flux ratios. However, 
anomalous flux ratios can also be caused by stellar microlensing if the source is small enough. One solution 
is to observe the narrow line region emission lines, where the lensed quasar emission is sufficiently extended 
to be unaffected by microlensing. Large samples of thousands of lenses could achieve sensitivity down to 
10® — 10^ solar masses required to probe the nature of Dark Matter. Unfortunately, at the moment only 3 
dozen quads are known, and only a handful of those are bright enough to be observable from the ground 
(e.g. Chiba et al 2005). 

Discovering thousands of quads will require surveys covering thousands of square degrees at sub- 
arcsecond resolution (Oguri & Marshall 2010). Between them, Euclid, WFIRST and LSST should discover 
over 10,000 lensed AGN, one sixth of which will be quad systems. Euclid and WEIRST have the potential to 
observe of order a thousand of these quads spectroscopically themselves, using narrow near infrared emis¬ 
sion lines to avoid the microlensing and enabling the wholesale application of the lens flux ratio experiment 
(Nierenberg et al. 2014). 

4.3 Deep drilling fields 

The ESST collaboration plans several “deep drilling” fields. While this aspect of the ESST survey is still 
under design, the broad plan is to reach AB ss 28 in ugriz bands (and shallower in y), which is comparable 
to the medium depth planned for the WEIRST SN survey (J=27.6 and H=28.1 over 8.96 deg^). The WFIRST 
deep SN fields will be deeper with depths of J=29.3 and H=29.4 over 5.04 deg^. ESST has currently 
identified four selected fields, each with areas of 9.6 deg^: EEAIS SI, XMM-ESS, COSMOS, and Extended 
Chandra Deep Field South. With three days of observation per band, WFIRST can cover each of these deep 
fields to AB = 28. 

The goals of these fields (Gawiser et al., white paper; Ferguson et al., white paper) are to (1) Test and 
improve photometric redshifts critical for ESST Main Survey science. (2) Determine the flux distribution 
of galaxy populations dimmer than the Main Survey limit that contribute to clustering signals in the Main 
Survey due to lensing magnification. (3) Measure clustering for samples of galaxies and Active Galactic Nu¬ 
clei (AGN) too faint to be detected in the Main Survey (4) Characterize ultra-faint supernova host galaxies. 
(5) Characterize variability-selected AGN host galaxies. (6) Identify of nearby isolated low-redshift dwarf 
galaxies via surface-brightness fluctuations. (7) Characterize of low-surface-brightness extended features 
around both nearby and distant galaxies. (8)Provide deep “training sets”, for characterizing completeness 
and bias of various types of galaxy measurements in the wide survey (e.g. photometry, morphology, stellar 
populations, photometric redshifts). All of these goals will be clearly enhanced by complementary data in 
the infrared. As already outlined in the report, infrared data will significantly improve the certainty of photo- 
z’s. WFIRST high spatial resolution will be particularly powerful in these deeper fields where confusion is 
likely to be significant. The deblending, color correction and morphology arguments in sec 3 apply even 
more strongly at the greater depth reached in these fields. 

Indeed, WFIRST provides great synergy by providing a high-resolution detection image that will allow 
ESST to go past the ground-based confusion limit. This in turn will motivate ESST to go deeper in gri than 
currently planned for the Deep Drilling Fields; by matching the ESST Main Survey filter time distribution, 
ESST will reach deeper than AB~29 in these three filters rather than the current baseline of AB~28.5, 
providing a better match to the deepest WFIRST imaging. 

It is worth noting that the VISTA VIDEO deep NIR survey regions are located at the centers of the 
already-approved ESST DDEs, and these will generate initial spectroscopy and multi-wavelength coverage 
even before WFIRST launches, making them optimal regions for covering deeper in NIR. However, if 


22 



WFIRST is able to cover all (or most) of the 9.6 deg^ LSST field-of-view, LSST can respond by skipping 
the significant fraction of time spent in y band in the current LSST ’’near-uniform DDF” approach. Using 
WFIRST imaging in lieu of LSST imaging in y band (and possibly z band) also offers a big improvement 
in both resolution and depth, as LSST only reaches AB=28.0,27.0 in z and y in the ’’near-uniform DDF” 
approach spending 60% of DDF observing time in those filters; significantly deeper ugri images could be 
obtained if this time is saved. Full areal coverage by WFIRST of as many LSST DDFs as possible, in 
zyJH, will motivate LSST to devote its observing time in those fields fo fhe deepesf possible ugri imaging, 
enabling a wide range of science wifh fhese ”ullra-deep drilling” fields. 

4.4 Coordinated supernovae observations 

If fhe WFIRST SN medium and deep surveys largel fhe same fields as one of fhe LSST deep drilling fields, 
fhen coordinated observafions would yield 8 band supernova lighf curves. The WFIRST IFU measuremenfs 
of fhese supernovae would yield a powerful fraining sef fhaf could greafly enhance fhe value of fhe large 
LSST SN sample. These coordinated observations would nol require any increase in observing lime by 
eilher projecl bul would enhance fhe refurn of bofh sefs of observations. 

4.5 Variable universe: follow-up on novel targets 

LSST’s cadence and depfh will allow if fo discover new classes of variable objecfs and “golden” fransienfs, 
examples of known classes of objecfs lhal have nearly ideal properfies for particular aspecfs of astronomical 
study. WFIRST’s resolution, spectral coverage, and wavelength coverage will make it a premier facility for 
following up on some of these objects in the near infrared. Since studies of individual targets do not take 
advantage of WFIRST’s large field of view, large ground-based telescopes wifh high qualify adaptive opfics 
will also be imporfanf fools for follow-up observafion. 

There are a hosf of pofenfial fargefs for follow-ups. Examples include: 

• “Macronovae” produced by fhe coalescence of binary neufron sfars (Li & Paczynski 1998, Kulkarni 
2005, Bildsfen el al. 2007). These explosive evenfs likely peak in fhe near-infrared and would be 
pofenfial fargefs for bofh WFIRST imaging and IFU follow-up. 

• LSST and WFIRST will be our mosf powerful combination of felescopes for defecfing afterglows of 
gravifafional mergers (Kasliwal and Nissanke 2014; Gehrels el al. 2014). Given fhe LIGO’s large 
error box for each evenl, Ihere will be many variable objecfs found af fhe LSSTAVFIRST deplhs, a 
program of join! moniloring (perhaps, combined wifh fhe deep field program described above) will 
enable a heller characlerizalion of fhe rich zoo of fransienfs. 

• Gamma Ray Bursls (Rossi el al. 2007). 
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5 Data management and analysis challenges 


The sections above make clear that in essentially every area of astronomy and cosmology, joint analysis 
of the LSST, WFIRST and Euclid data will provide significant benefits. These will be made possible by 
linking the data from the different surveys and providing a common access point for interrogating the data in 
a user-friendly way with the appropriate tools. This would enable scientists to explore connections between 
the data sets and stimulate novel community research leveraging the combined data. The publication of 
detailed, accurate documentation along with the cross-referencing and associated data releases would be 
the most effective way to encourage and support the community in pursuing such research. This section 
describes the issues involved in carrying out such a program for the photometry and simulations required to 
analyze and interpret the data. 

5.1 Multi-resolution photometry 

There are quite a few challenges to producing accurate, consistent photometry across different filters with 
widely different resolutions, as would be the situation for combining LSST data with either WFIRST or 
Euclid. 

We can assume that we will have a catalogue of the positions of objects detected in at least one band 
(including single-band detections in u or H), and that we know the PSF in each band. There are issues with 
respect to the chromatic effects within a single band as well as spatial variation of the PSF across the field of 
view, which are imporfanf, buf we can neglecf fhem for fhis discussion. We also assume fhaf blending issues 
have been faken info accounf, so we have a sef of pixels in each band fhaf correspond fo a single objecf. 

For poinf sources, phofomefry is relafively slraighfforward. Because all poinf sources have an idenfical 
profile (again, ignoring for now fhe spatial variation of fhe PSF) we may use any aperfure fo measure fhe 
relative flux of all poinf sources in an image (a PSF aperfure is opfimal). We may fhen esfimafe fhe image’s 
zero-poinf which gives us fofal fluxes for all sources. We may repeal fhis in all fhe bands, resulfing in a 
cafalogue of fluxes and colors. 

Measuring fhe fluxes of galaxies is harder. If is nol possible fo simply measure fhe fofal flux in an 
aperfure large enough fo include fhe enfire galaxy as fhe measuremenfs are foo noisy, and fhis is nol even 
a well-defined concepl. On fhe olher hand, a small aperfure measures a differenl fraction of fhe flux of 
differenl galaxies, so we cannol simply apply an aperfure correclion as in fhe slellar case. Even if a galaxy 
has fhe same profile in each band Ihese fractions are differenl if fhe PSF is band-dependenl; if fhe profiles 
vary fhe difficulties are multiplied. 

If we can model galaxies adequalely, fhis alleviates fhe difficullies of flux measuremenl. The PSF- 
convolved model flux is a reasonable definition of fhe galaxy’s fofal flux in each band. However using 
model-based phofomefry is nol a panacea, as if fhe model musl be correcl. For example, choosing a model 
in one band and applying if in anolher will nol yield fhe correcl color for a galaxy wilh a red bulge and a 
blue disk, a scale lenglh lhal’s a funcfion of band, spiral arms, HII regions, or any olher kind of subslruclure 
wilh a differenl color from fhe overall galaxy lighl. 

If fhe seeing is (or can be made) fhe same in all bands we are allowed fo choose a single model. The re¬ 
sulting fluxes are nol fhe galaxy’s fofal fluxes, as differenl componenfs are nol necessarily weighfed correcfly, 
buf Ihese “consislenf fluxes” do represenl fhe fluxes of a cerlain sample of fhe galaxy’s slellar populalion. 
However, degrading all bands fo a lowesl common denominalor would deslroy much of fhe exlra informa¬ 
tion presenl in fhe higher-resolulion space dala. Determining fhe opfimal way fo esfimafe consislenf fluxes 
white keeping as much of fhis high-resolulion information as possible is still an open research question. 

Given fhe complexify of fhe problem of measuring galaxy phofomefry across a decade in wavelenglh 
wilh varying PSFs, Ihere is no opfimal approach. Differenl science questions will require differenl ways of 
combining fhe mulli-color dala. 
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• Photometric Redshifts: We may use an objeet’s total or eonsistent fluxes to find a loeation in multi- 
eolor spaee and henee its redshift probability funetion, p(z). If p{z) is unimodal and sharp, the two 
estimators prediet nearly the same photo-z. However, in general this will not be so and p{z) will be 
different for the two ehoiees of flux. 

• Choice of Galaxy Model: In the rest-frame optieal, a simple Sersie or eonstrained bulge-disk model 
is probably suffieient to return reliable fluxes (sueh models down-weight the eontribution from lo¬ 
calized star-forming regions, changing p(z) at least in principle). However, modeling galaxies across 
the broad redshift range covered by WFIRST, Euclid, and LSST will likely require more complex 
models (including, for example, the uv-bright knots that contribute the bulk of the flux in the far-uv). 
Constraining these models is going to be hard; one approach would be to demand that all components 
have a consistent p{z). 

• Photometry of Blended Objects If our models indeed accurately describe our galaxies then the de¬ 
blending problem is no different from the galaxy photometry problem; instead of fitting multiple 
components to a single galaxy we fit multiple components to multiple galaxies, relaxing but not aban¬ 
doning the constraints on components’ p{z). In practice it is not clear how well this will work, and 
there are intermediate schemes (similar to the deblender used in SDSS) that fit simplified models to 
the ensemble of objects to assign the flux to the child objects and then proceed with photometry. This 
is an unsolved problem, but given the upcoming data, we will need to must significant progress in the 
coming years. 

5.2 Joint simulations 

WFIRST, Euclid and ESST will generate a remarkable set of observations in the optical and near-infrared 
wavebands. Their cosmological interpretation, especially in combination, will be very powerful but also 
very challenging. Many systematic effects will have to be disentangled from fundamental physics to fully 
exploit the power of the measurements. Simulations of all aspects of the experiment - at the cosmological 
scale, of the instruments, and of the data processing and analysis software - are critical elements of the sys- 
tematics mitigation program. Joint analyses of multiple surveys require these simulations to have consistent 
interfaces, to enable the same realizations of cosmological volumes or multi-wavelength simulated galaxies 
to be passed through instrument-specific simulators, and ultimately to ensure that the correct correlations 
and joint statistical properties are represented when mock analyses are performed on the simulated output 
catalogs. 

We focus here on the role of coordinated simulations in two areas - cosmological simulations and in¬ 
strument/pipeline simulations. 

5,2.1 Cosmological simulations 

Cosmological simulations play key roles within the interpretation of large scale structure surveys: (i) they 
provide controlled testbeds for analysis pipelines and tools for survey design and optimization via sophis¬ 
ticated synthetic skies populated with galaxies targeted by the specific survey, (ii) they provide predictions 
for different fundamental physics effects, including dynamical dark energy, modified gravity, neutrinos, and 
non-trivial dark matter models, (iii) they provide modeling approaches for astrophysical systematics, includ¬ 
ing gas dynamics, star formation, and feedback effects, and (iv) they provide important information about 
the error estimates via covariance studies. We will briefly elaborate on the synergies between the simulation 
programs required for WFIRST and ESST in these four areas in the following. 

Synthetic sky maps — In order to build synthetic sky maps for Euclid, WFIRST and ESST, large-volume 
A^-body simulation with very high force and mass resolution are essential. Next, the simulations have to be 
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populated with galaxies using sub-halo abundanee matehing or semi-analytie modeling teehniques, tuned 
to mateh the observed population of galaxies. While the wavelength ranges of Euelid, WFIRST and LSST 
will be different, a eommon eode infrastrueture will be benefieial in some areas (e.g., in exploring model 
parameter spaee) and eritieal in others (e.g., simulating realistie merged multi wavelength eatalogs). 

Precision predictions for fundamental cosmological statistics — The preeision requirements for eos- 
mologieal observables and the fundamental physies to be explored is very similar for Euelid, WEIRST and 
ESST and a joint program in this area would be very fruitful. Both surveys are promising observations 
of large seale strueture measurements at the sub-pereent level, a tremendous ehallenge for eosmologieal 
simulations. Besides the preeision ehallenge, different fundamental physies effeets have to be explored 
systematieally for both surveys. 

Astrophysical systematics — Some of the systematies that have to be aeeounted for to interpret the 
observations from Euelid, WEIRST and ESST are eommon, sueh as the effeets of baryonie physies on the 
weak lensing shear power speetrum. Sinee the aeeuraey requirements are similar, a joint program would be 
benefieial to both surveys. However, there are signifieant differenees beeause of the differing ground and 
spaee-based nature of the two programs, different instrumentation and assoeiated wavebands, and different 
galaxy populations and densities. Nevertheless, the underlying sky eatalogs ean be generated from the same 
set of simulations. 

Covariance matrices — Sinee the Euelid, WEIRST and ESST surveys will have overlapping volumes, 
eombining eosmologieal eonstraints from their data sets will require understanding the full eovarianee ma¬ 
trix of all observables extraeted from the data sets, using simulations of the entire relevant volume. While 
the aeeuraey requirements for these simulations are mueh less stringent than for other applieations, the large 
number of simulations required will represent a eomputational ehallenge. 

5.2.2 Instrument and pipeline simulations 

Instrument and pipeline simulations are also a eritieal element of any preeision large seale strueture program: 
they inform the relationship between the astronomieal “seene” in the Universe and the eatalogs, parameters, 
and other information that is reeorded in the proeessed data. They also provide an opportunity to exereise the 
data reduetion and analysis tools prior to first light. The instrument and pipeline tools must be eonsidered 
together, sinee in a preeision experiment the pipeline is an integral part of the measurement proeess. Unlike 
eosmologieal simulations, these tools are highly instrument-speeifie: for example, the atmosphere is only 
ineluded for ground-based observations, the patterns of ghosts and diffraetion features are eompletely dif¬ 
ferent for the ESST, Euelid and WEIRST eonfigurations, and internal optieal effeets, eross-talk, and readout 
artifaets are fundamentally different in silieon CCDs versus NIR arrays. 

Despite the differenees in instrument simulators, it is eritieal that they be able to take eompatible input 
data. Preeision measurements of galaxy elustering and weak lensing depend on not just eounting objeets and 
measuring shapes, but understanding seleetion biases, blending effeets, and error distributions of measured 
parameters. The eombination of WFIRST, Euelid and ESST data will be used to measure photometrie 
redshifts, and for both objeet-by-objeet and statistieal eomparisons of number density and shear. Simulation 
of these applieations requires that the same objeets be fed through both pipelines, and henee that the input 
seene (simulated stars and galaxies) be eonsistent and have the proper eorrelations of objeet properties aeross 
wavebands. While proprietary pipelines do not need to be shared outside of the individual eonstortia, the 
inputs of these pipelines must remain eompatible. 
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6 Conclusion 


The scientific opportunity offered by the combination of data from LSST, WFIRST and Euclid goes well 
beyond the science enabled by any one of the data sets alone. The range in wavelength, angular resolution 
and redshift coverage that these missions jointly span is remarkable. With major investments in LSST and 
WFIRST, and partnership with ESA in Euclid, the US has an outstanding scientific opportunity to carry out 
a combined analysis of these data sets. It is imperative for us to seize it and, together with our European 
colleagues, prepare for the defining cosmological pursuit of the 21st century. 

The main argument for conducting a single, high-quality reference co-analysis exercise and carefully 
documenting the results is the complexity and subtlety of systematics that define this co-analysis. Falling 
back on many small efforts by different teams in selected fields and for narrow goals will be inefficient, 
leading to significant duplication of effort. 

For much of the science, we will need to combine the photometry across multiple wavelengths with 
varying spectral and spatial resolution - a technical challenge. As described in Section 3, the joint analysis 
can be carried out in ways that have different computational demands. The most technically demanding joint 
analysis is to work with pixel level data of the entire area of overlap between the surveys. Many of the goals 
of a joint analysis require such a pixel-level analysis. If pixel-level joint analysis is not feasible, catalog-level 
analysis can still be beneficial, say to obtain calibrations of the lensing shear or the redshift distribution of 
galaxies. Hybrid efforts are also potentially useful, for example using catalog level information from space 
for deblending LSST galaxies, or using only a mutually agreed subset of the data for calibration purposes. 
However the full benefits of jointly analyzing any two of the surveys can be reaped only through pixel-level 
analysis. 

The resources required to achieve this additional science are outside of what is currently budgeted for 
LSST by NSF and DOE, and for WFIRST (or Euclid) by NASA. Funding for this science would most 
naturally emerge from coordination among all agencies involved, and would be closely orchestrated scien¬ 
tifically and programmatically to optimize science returns. A possible model would be to identify members 
of the science teams of each project who would work together on the joint analysis. The analysis team would 
ideally be coupled with an experienced science center acting as a focal point for the implementation, and 
simultaneously preparing the public release and documentation for broadest access by the community. 
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